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Abstract 

Baryons containing two heavy quarks are treated in the Born-Oppenheimer ap- 
proximation. Schrodinger equation for two center Coulomb plus harmonic oscillator 
potential is solved by the method of ethalon equation at large intercenter separa- 
tions. Asymptotical expansions for energy term and wave function are obtained 
in the analytical form. Using those formulas, the energy spectra of doubly heavy 
baryons with various quark compositions are calculated analytically. 
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Introduction 



The investigation of properties of hadrons containing one or more heavy 
quarks is very important for understanding the dynamics of quark and 
gluon interactions. Presently at the LHC, B-factories and the Tevatron 
with high luminosity, several experiments have been proposed, in which a 
detailed study of baryons containing two heavy quarks can be performed. 
In particular, in the forthcoming experiment at CERN, the COMPASS 
group is going to find the doubly charmed baryons and study their phys- 
ical properties. In this connection, doubly heavy quark baryons are now 
becoming one of the most exciting subjects in particle physics. Therefore, 
theoretical predictions of properties of doubly heavy quark baryons acquire 
a big significance for the forthcoming experimental study of these particles. 
So far there have been various approaches by which their mass spectra and 
other properties can be calculated. One of them is the nonrelativistic quark 
model which gives relatively accurate results for baryon spectra Ql], @, |3[]. 
The possible quark compositions of doubly heavy quark baryons are ccg, 
cbq and bbq, where q denotes a light w, d or s quark. Note that the baryons 
containing the top quark(s) are not practical subject here because a top 
quark is extremely heavy and hence we have no chance to find them as 
stable hadrons. The doubly heavy quark baryons may be considered as 
an analogue of the hydrogen molecular ion ', which has been treated 
successfully in the Born-Oppenheimer approximation. The same approx- 
imation is expected to be efficient even for doubly heavy quark baryons, 
though there exist some differences between these baryons and systems. 
One of them is, for this case, the appearance of the confining potential in 
addition to the QCD Coulomb potential. As is well known, the variables of 
the Schrodinger equation with two-center Coulomb plus confining poten- 
tial cannot be separated for their kinematical variables, in general. To our 
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knowledge, the two-center potential which allows the separation of vari- 
ables is only the two-center Coulomb plus harmonic oscillator potential. 

In this paper, we treat QQq baryons in the nonrelativistic approach by 
using the solution of the Schrodinger equation with two-center Coulomb 
plus harmonic oscillator potential, i. e. the well-known method of ethalon 
equation which is widely used for solving Schrodinger equation with two- 
center pure Coulomb potential in the physics of H 2 [H, 0, |J. First, we 
give a general scheme of treatment QQq baryon in the Born-Oppenheimer 
approximation. Then, the two-center Schrodinger equation with two-center 
Coulomb plus harmonic oscillator potential is analytically solved with some 
approximation: the energy term of the light quark moving in the field of 
two heavy quarks is obtained in the form of asymptotical expansion over 
the inverse power of the distance between heavy quarks. Finally, we give 
an analytical formula of the baryon energy spectrum for QQq. 

Doubly heavy quark baryon in the Born-Oppenheimer 
approximation 

In the Born-Oppenheimer approximation the wave function is split into 
heavy- and light-quark degrees of freedom 

«r) = Ei(%Wr), 

n 

where R is the distance between two heavy quarks and r is the distance 
between light quark and center-of-mass of the heavy-quark pair. The light 
quark wave function ^(r, R) and its energy term E(R) can be found from 
the Schrodinger equation 

[-^_A + V(n) + V(r 2 )]i/j = E(R)i/j, 

where r\ and r 2 are the distances between light and heavy quarks, Q\ and 
Q2, respectively. The binding energy of this system is approximated by 
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the equation: 

1 -A + V QQ (R)+E(R)]<f> = e<f> 1 



2Mg Q 

where Mqq is the reduced mass of QQ. 

A quark potential with Coulomb plus harmonic confinement for this 
baryon is given by[ 



Vinj) = X)y \\MVo - Ar% + Hi) = -^(Vo - Ar% + 

In the field of two heavy quarks with this potential, the motion of a light 
quark can be nonrelativistically described by the following Schrodinger 
equation, 

[-1a - - - - + u\r\ + r\) - \v<M> = E{R)^ (1) 

where Z = 2a s /3 and uj 2 = 2A/3. 

In the prolate spheroidal coordinates defined as 

Z = (K e < oo), rj = (-Krj< 1), 

the potential term in eq. ([!]) can be written in the form 

2 q(fl+6fa) q^tf 4 

%r2) -"jpT^r + i-^ ( 2 ) 

where 

a(0 = 2ZR - ^fea 2 - 1), b( V ) = 2ZR - - 1). 

As is well known||], the Schrodinger equation with the potential in the 
form of eq.(0) is separable in the prolate spheroidal coordinates. Then it 
is convenient to use 

U{0 V{rj) e ±im<t> 

where and m are azumuthal angle and azimuthal quantum number, re- 
spectively. After substituting this into eq. ([!]), we obtain from the following 
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ordinary differential equations connected with separation constants A and 
rn: 



u"(() + £ + h -f^ - r^e + ^Mf) = o, (3) 



.h 2 hX ,42 1 — m 2 



where a = 2Z/\ / 2E' and 7 = cu 2 /8E /2 , and further 



/i = v^fl, (5) 

with 

Finiteness and continuity of the wave function ^ in the whole space lead 
to the following boundary conditions for the functions U and V: 

U(0k=i=0, 17(0 le— k» — ► 0, (6) 

V(v) \ v =±i= 0. (7) 



Asymptotics of quasi-angular equation 

We will approximately solve eqs.(|3|) and(g) for large R by the method of 
ethalon equation. This method is successfully applied to the solution of 
nonrelativistic two center Coulomb problem [||, |5|, |5| and in the theory of 
diffraction of waves. Details on the method of ethalon equation are given 
in ||, §, §, 0| , also briefly described in Appendix. 

Let us start from the angular equation of (§j). As an ethalon equation 
for eq.(f§), we choose the following Whittaker equation [|8j: 

TT w/ r h A h 2 k 1 — m 2 lTTr 

W" + — — + + —r^-]W = 8 

4 z 4z 2 
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and seek a solution in the form 

V=[z'(rj)]^M Kf (h 2 z), (9) 



where Mk^(h 2 z) is the solution (regular at zero) of eq.([8|). Substituting 
(0) into flD and taking into account (§), we get the following equation for 
z: 

z' 2 . 1N2 1,1 kz' 2 A 
T- 7(X - 1} -7^ ( 4 + ^-2x(l-x/2) ) + 

^^i^y-^-i^"^ ' (10) 

where 

3 z" z w 

and t = J" , x = 1 + 77. 

Requirement of coincidence at the transition points [0, 0] 

z(z) U=o= 0, 

leads to the following "quantum condition" 

A ^ 2^(0) + f - 1]. (11) 

We will seek the solution of eq.( |TD|) and eigenvalues A in the form of fol- 
lowing asymptotical expansion: 

_ °° Zk_ ~ Xk 

z ~ L^i 7 * ~ L^i 7 • 

k=0 n fc=0 a 

Substitution of these expansions into (jXDj) gives us the recurrence system 
of differentil equations for z : 



4 = 2 7 l(x-l), 



z[ = 0, 

z , _ 1 2fc2j (^) 2 2A _ z£ 

" 2 24 z 24 4x(l-x/2) 2 
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and for A : 

A = 2*4(0), 
Ai = 2kz[(0), 

Aa = 2*4(0) + mS"1), 



Solving these recurrence equations, we obtain 

A<"» = 4^ + ^^ + 0(1), (12) 

for A, and 

2 = 1 \ x {2 - x ) + -roPln{l -x) + 0(-?-), (13) 

for z. 

^From boundary conditions one can obtain for quantum number 0: 

m + 1 

where g = 0, 1, 2, 

Asymptotics of quasi-radial equation 

As an ethalon equation for eq.(^), we take the following equation : 

W" + [h 2 s - h A y 2 - = 0, (14) 

a solution of which is expressed by the confluent hypergeometric functions [§], 

a 

c _^V EV S-2C-1 1 4 2 



where c = 1+ ^ 2 + 3 , 
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Boundary condition (|3j) and properties of functions F || give rise to the 
following expression for s: 

s = 4n + Vm 2 + 3 + 2. 

Substituting 

U=[y(0]^W(y(0) 

into eq.(|3|), we obtain 

y 2 y' 2 2 l_ 1_ /2 A 1 o£ 4r 3-4r ?/ 2 1_ 

(15) 

After substitution of 

= ^, 

^ 4 ' 

this equation can be reduced to the form 

,/a , \ 1,1, 1 N /2 A . 1 a(* + l) r ,0 /2 4 , r , n „ 

(16) 

where £ = £ — 1. The quantization condition which follows from 4>{x) = 
is written in the form 

A = -2^(0) + I - i[ J + 1] |« . (17) 
Inserting the asymptotical expansions 

~ fa oo Afc 

fc=0 a fc=0 " 

into eq. fllfip and solving the equations obtained herewith, we get the fol- 
lowing result 

y = 2 7 i(t 2 + 2tf* + tUt"^ 2 + 2t)~hn(t + 1) + 

at 

la 7 -f(t 2 + 2«)-^n^>+0(i), (18) 



for y, and 



for A. 



. (t\ i a 4r — s£ scry 4 ^, 1 < , * 

x( ° = -h + — -^r + (19) 
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Asymptotical expansion for energy 



Asymptotical expansions ( |I2|) and ([T9]) give us an expression for the energy 



term in the form of multipole expansion. In order to obtain this expansion 
one should insert 

*-* + ! + ! + ... 

into eqs.( |I2|) and (|I9j). Equating to be A^ and taking into account 



we get the following equations for coeficients Ei, E<i, 



Ei = ^[(su - 2ku- 1 )(2E )-2 + (4s 2 - 16k 2 - 16r)(2 J B ) i ], 
bZ 

E 2 = -El + 2suj- 1 Eq + E 1 (2E )5Z" 1 (16r 2 + 16k 2 - 4s 2 ), 
2 



Now we need to find Eq. In order to find this value, we note that for 
R — ► oo, E' = Eq and hence we have 

E = E Q + ^- + -V . (20) 
On the other hand, for large R we have 

97 oo r E>2 

V(n, r 2 ) = £(-=)' l Pi{cos9) + to 2 [{r 2 + 2rRcos6 + — ) + 
R i=o R 4 

(r 2 - 2rRcos6 + ^-)] » u; 2 (2r 2 + - jV . (21) 

Hence, for the energy term with this potential we obtain 

E = MN + \) + ^- + \v , (22) 

where N = n J rq + m+ l is the principial quantum number. Comparing 
eqs.(gy) and (|22|), we obtain 

E = 2u(N + ^). 
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Thus, the following asymptotical expansion is obtained for the energy 
term of light quark in the field of two heavy quarks: 

4 cu 2 R 2 „ Ei E 2 
E = - 3 V " + — + E " + R + IV + - 

QQq baryon spectra 

As mentioned above, the QQq binding energy can be finally obtained by 
solving the Schrodinger equation 

[-ttt^A + V QQ (R) + E(R)]<f> = ecji. (23) 

If one takes E(R) in the form 

E = — Vo H h E H , 

3 2 K 



for 



V QQ (R) = cu 2 R 2 - | - ?^ , 



then eq. (fZ3f) can be rewritten as 



where Z' = Z - E u J 2 = \u 2 and Vj = 2V - #o- 

To solve this equation, we use the result of |TT| where a method for an 



analytical solution of the Schrodinger equation with potential 

V(R) = -^ + \R k 
R 

was offered. Details on this method and its application to our potential 
are given in Appendix. Application of this method to eq. ( gjp gives us 



e Nnl = 2u{N + |) + [Z' 2 u' 6 r nl ] ^ - 2V , (25) 

where N is the principial quantum number of the light quark moving in 
the field of QQ, r n \ is defined in Appendix. The formula ( fZ5|) describes the 
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energy spectrum of the QQq baryon. In tables 1, 2 and 3, the mass spectra 
of ccg, bbq and bcq baryons calculated using the formula fl25Q are given, 
respectively. The following values of potential parameters are choosen in 
this calculation a s = 0.39, to 2 = 0.174GW 3 , V = O.ObGeV for the 
potential 



Conclusion 

In this work we have treated doubly heavy baryons in the Born-Oppenheimer 
approximation. The following two problems have been solved in the frame- 
work of this approximation: (l)the Schrodinger equation for two-center 
Coulomb plus harmonic oscillator potential and (2)the Schrodinger equa- 
tion for central symmetric Coulomb plus harmonic oscillator potential. As 
the final result an analytical formula for the energy spectrum of baryons 
containing two heavy quarks is derived. Obtained formula is applied for 
the calculation of mass spectra of doubly heavy quark baryons with vari- 
ous quark compositions. The above analytical results could be useful for 
further numerical calculations in non-asymptotical region. 

Appendix 

1 The scaling variational method and its application 
to Coulomb plus confining potential 

Consider the following Hamiltonian 




H = — A + V(r) 



(26) 



which obeys the eigenvalue equation 



Hljj n i = Enil/Jni, < Iprd \ ^n'V >= <W$F 



(27) 
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n,n' = l,2,.., 1,1' = 1,2,.., 

where n and / denote the principal and angular quantum numbers, respec- 
tively. To solve this equation, we start from a set of functions {4> n i} which 
are the eigenfunctions of an arbitrary central field Hamiltonian Hq: 

Ho4> n i = e n i(j) n i, < 4> n i | (f) n >i> >= 5 nn r5u>. (28) 

Then, we construct the functionals 

e nl {a) =< 05 | Hfa >, (29) 

where 

<\> a nl = a i/2 4) nl (ar). 
The value for the a is determined by 

(^)(a = a)=0. (30) 

Let us now apply this method to our problem, i.e. to the Schrodinger 
equation with potential 

Z' 

V(r) = -— + u' 2 r 2 . 
R 

As Hq, we choose pure Coulomb Hamiltonian, i.e. 

H = -Ia - Z -. 
u 2 R 

Then, e n i = —§^z with n = n r + / + 1, where n r is the radial quantum 
number. 

According to the above procedure, we have 



s nl (a) =< C | H | C > + < C | | ^ >= (31) 

+ [5n 2 + l-3/(/ + l)], 



2n 2 
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For calculation of the second matrix element, we have used the well-known 
expression for average value f 2 in Coulomb field which is given in |T1[ . 
^From de(a)/da = 0, we obtain 

9, ) 2 <n 4 

ao = { -^- [5n2 + 1 - 3 ' (i+1)1}1/5 ' 

Then, for the energy level one can get the following analytical formula 

= [ZVrn,] 1 / 5 , (32) 

where 

r^ = n 2 [5n 2 +l-3/(/ + l)] 3 . 
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2 The formal procedure of the method of ethalon 
equation 

Let's consider the following second order differential equation: 

y"(x)+p 2 [\-q(x)}y(x) = (33) 

in the interval [a, b] . Let in this interval eq. ( |33| ) has one transition point (poles 
and zeros of function Q(x, A) = q(x) — A are called transition points of this 
equation) . 
Equation 

w"(z) -p 2 R(z)w(z) = (34) 

which has the same or close transition points as eq.([J4]) is called the ethalon 
equation for eq . (|33|) . 



Solution of eq.(|33|) we will seek in the form 

y(x) = [z'(x, p)]~ 1 / 2 w(z(x,p)) (35) 

where w is the solution of eq . (|34Q . Inserting ( |35|) in to eq. fl33Q and tak- 
ing into account eq.(^) we obtain the following (nonlinear) differential 
equation for z(x,p): 

R{z)z' 2 - Q(x, A) - x} = (36) 

2p z 

here 

3 z" z'" 

{*,*} = -2 { v f+ -' 

In the case of eq. (f|) we have for Q 



_ r h 2 hX 7 4 2 1 - m 2 
and for R (from eq.(g) 

r h A h 2 k 1-m 2 - 

R = — — T + + A 9 

4 z 4z 2 
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So, for z one obtains eg . (jTDj) . 
T -7(-D 2 -^(- 



A 



)+ 



T 



Z 
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z 2x{l-x/2) J h? K x 2 {l 



2h 2 



{z, x} = 
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Table 1. The mass spectrum of ccq baryon (in GeV) calculated using 
the formula fl25p ; m q = 0.385 GeV, m c = 1.486 GeV, n\ and n d are the 
principial quantum numbers of light quark and cc diquark, respectively, L 
is the orbital quantum number of cc diquark. 





ni,n d ,L 


mass 


m,n d ,L 


mass 


ni,n d ,L 


mass 


ni,n d ,L 


mass 




1,1,0 


3.661 


1,1,1 


3.613 


1,2,2 


3.649 


1,3,3 


3.694 




1,2,0 


3.730 


1,2,1 


3.708 


1,3,2 


3.764 


1,4,3 


3.825 




1,3,0 


3.816 


1,3,1 


3.799 


1,4,2 


3.872 


1,5,3 


3.949 




1,4,0 


3.914 


1,4,1 


3.901 


1,5,2 


3.988 


1,6,3 


4.077 




1,5,0 


4.024 


1,5,1 


4.012 


1,6,2 


4.110 


1,7,3 


4.210 




2,1,0 


3.839 


2,1,1 


3.791 


2,2,2 


3.828 


2,3,3 


3.873 




2,2,0 


3.908 


2,2,1 


3.887 


2,3,2 


3.942 


2,4,3 


4.003 




2,3,0 


3.994 


2,3,1 


3.978 


2,4,2 


4.051 


2,5,3 


4.127 




2,4,0 


4.093 


2,4,1 


4.079 


2,5,2 


4.166 


2,6,3 


4.255 




2,5,0 


4.202 


2,5,1 


4.190 


2,6,2 


4.289 


2,7,3 


4.388 




3,1,0 


4.018 


3,1,1 


3.970 


3,2,2 


4.006 


3,3,3 


4.051 




3,2,0 


4.086 


3,2,1 


4.065 


3,3,2 


4.120 


3,4,3 


4.182 




3,3,0 


4.172 


3,3,1 


4.156 


3,4,2 


4.229 


3,5,3 


4.305 




3,4,0 


4.271 


3,4,1 


4.257 


3,5,2 


4.344 


3,6,3 


4.433 




3,5,0 


4.380 


3,5,1 


4.369 


3,6,2 


4.467 


3,7,3 


4.567 
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Table 2. The mass spectrum of bbq baryon (in GeV) calculated using 
the formula fl25Q ; m q = 0.385 GeV, = 4.88 GeV, n\ and n d principial 
quantum numbers of light quark and bb diquark, respectively, L is the 
orbital quantum number of bb diquark. 



ni,n d ,L 


mass 


ni,n d ,L 


mass 


ni,n d ,L 


mass 


ni,n d ,L 


mass 


1,1,0 


9.890 


1,1,1 


9.874 


1,2,2 


9.886 


1,3,3 


9.900 


1,2,0 


9.911 


1,2,1 


9.905 


1,3,2 


9.922 


1,4,3 


9.942 


1,3,0 


9.939 


1,3,1 


9.934 


1,4,2 


9.957 


1,5,3 


9.981 


1,4,0 


9.970 


1,4,1 


9.966 


1,5,2 


9.993 


1,6,3 


10.022 


1,5,0 


10.005 


1,5,1 


10.001 


1,6,2 


10.032 


1,7,3 


10.064 


2,1,0 


10.096 


2,1,1 


10.081 


2,2,2 


10.093 


2,3,3 


10.107 


2,2,0 


10.118 


2,2,1 


10.112 


2,3,2 


10.129 


2,4,3 


10.149 


2,3,0 


10.146 


2,3,1 


10.141 


2,4,2 


10.164 


2,5,3 


10.188 


2,4,0 


10.177 


2,4,1 


10.173 


2,5,2 


10.200 


2,6,3 


10.229 


2,5,0 


10.212 


2,5,1 


10.208 


2,6,2 


10.239 


2,7,3 


10.271 


3,1,0 


10.303 


3,1,1 


10.288 


3,2,2 


10.300 


3,3,3 


10.314 


3,2,0 


10.325 


3,2,1 


10.318 


3,3,2 


10.336 


3,4,3 


10.356 


3,3,0 


10.353 


3,3,1 


10.347 


3,4,2 


10.371 


3,5,3 


10.395 


3,4,0 


10.384 


3,4,1 


10.380 


3,5,2 


10.407 


3,6,3 


10.436 


3,5,0 


10.419 


3,5,1 


10.415 


3,6,2 


10.446 


3,7,3 


10.478 
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Table 3. The mass spectrum of bcq baryon (in GeV) calculated using the 
formula Q25|); m q = 0.385 GeV, m b = 4.88 GeV, m c = 1.486 GeV, m and 
n d principial quantum numbers of light quark and be diquark, respectively, 
L is the orbital quantum number of be diquark. 



ni,n d ,L 


mass 


ni,n d ,L 


mass 


ni,n d ,L 


mass 


ni,n d ,L 


mass 


1,2,0 


7.217 


1,2,1 


7.160 


1,3,2 


7.178 


1,4,3 


7.199 


1,3,0 


7.259 


1,3,1 


7.206 


1,4,2 


7.233 


1,5,3 


7.263 


1,4,0 


7.307 


1,4,1 


7.251 


1,5,2 


7.286 


1,6,3 


7.324 


1,5,0 


7.361 


1,5,1 


7.300 


1,6,2 


7.343 


1,7,3 


7.386 


2,1,0 


7.438 


2,1,1 


7.355 


2,2,2 


7.403 


2,3,3 


7.452 


2,2,0 


7.471 


2,2,1 


7.414 


2,3,2 


7.432 


2,4,3 


7.454 


2,3,0 


7.513 


2,3,1 


7.461 


2,4,2 


7.488 


2,5,3 


7.518 


2,4,0 


7.562 


2,4,1 


7.505 


2,5,2 


7.541 


2,6,3 


7.579 


2,5,0 


7.615 


2,5,1 


7.555 


2,6,2 


7.597 


2,7,3 


7.641 


3,1,0 


7.692 


3,1,1 


7.669 


3,2,2 


7.687 


3,3,3 


7.709 


3,2,0 


7.726 


3,2,1 


7.716 


3,3,2 


7.743 


3,4,3 


7.773 


3,3,0 


7.768 


3,3,1 


7.760 


3,4,2 


7.796 


3,5,3 


7.833 


3,4,0 


7.816 


3,4,1 


7.810 


3,5,2 


7.852 


3,6,3 


7.896 


3,5,0 


7.870 


3,5,1 


7.864 


3,6,2 


7.912 


3,7,3 


7.961 
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